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I review recent evaluations of the hadronic contributions to the muon anomalous 
magnetic moment (g — 2)fi and to the effective fine structure constant a{Mz)- A 
new estimate for the hadronic shift Aa'^^^^{M^) = 0.027782±0.000254 is presented 
which implies a~^(Mz) = 128.913 ± 0.035. It is based on a recent perturbative 
calculation of the Adler function which includes mass effects up to three-loops in a 
MOM scheme and requires little ad hoc assumptions. I then discuss perspectives 
for possible improvements in estimations of ajj'*'^ which we expect from the 
factory Daphne at Frascati. 



1 Status 

I briefly review the status of the estimations of hadronic vacuum polarization 
effects. I have presented a similar report some time ago in Ref. EJ and the 
present contribution should be considered as an Addendum to my previous 
summary. Vacuum polarization effects play a role in many places in physics 
(charge screening). Their precise knowledge is most important for the inter- 
pretation of high energy precision experiments where theoretical predictions 
depend on the effective fine structure constant a{s) = a/{l — Aa{s)) as well 
as for low energy precision experiments like for the anomalous magnetic mo- 
ment a^, one of the most precisely measured quantities in physics. The photon 
vacuum polarization amplitude Il'.^{q^) is defined by 

= * / < 0\Tj; {x) J2 (0) |0 > 

= n; (g") (1) 



"^Presented at the IVth International Symposium on Radiative Corrections, Barcelona, 
September 1998 (to appear in the proceedings, ed. J. Sola); 




and determines the full photon propagator 

i 1 , . 

as well as the shift in the fine structure constant 

Aa(g2) = e2(n;(0)-n;(g2)) . (3) 

The latter is subtracted at zero momentum, i.e., in the classical limit. At the 
one-loop order in perturbation theory one obtains 

with ^ the MS scale, Ncf the color factor, y — irnj/s and /3 — y/1 — y. 

The electromagnetic current is the sum of a leptonic and a hadronic part. 
The leptonic part can he calculated in perturbation theory and is known in 
QED up to three-loopsi: Aaicp = 0.031497687 (2-loop fraction - 2.5 x 10"'\ 
3-loop fraction - 3.4 x 10~^). The quark contribution cannot be calculated 
reliably in perturbation theory because of low energy strong interaction effects. 
Fortunately existing e~^e~ data allow us to evaluate the hadronic contributions 
by means of the dispersion integral 

where 

R.{s) - "^fr"^^*r^"f°"^^ = 12.1mn;(.) (6) 

has been measured in e~''e~-annihilation experiments up to about i?cut ~ 40 
GeV. By virtue of the asymptotic freedom of QCD the high energy tail can be 
safely calculated using pQCDB. 

Taking a conservative attitude we have compiled the existing e"'"e~-data (non- 
perturbative) and evaluated En. (||) for £ Ecut ^ 40GeV. After including 
the perturbative tail we foundQ 

Aa|f2i(M|) = 0.02804 ± 0.00064 (Eidelman, Jegerlehner 95) 
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Figure 1: Evaluations of l/a{Mz)- 



which leads to a^'^iMz) = 128.89 ± 0.09. It is obviously safe to use pQCD 
in the continuum above the ip resonances and below the T resonances (from 
5.0 GeV to Mr) and in the continuum above the T resonances above about 
12 Gevi. We then obtaini Aa'^^^iMj) = 0.02801 ± 0.00058. Utilizing the 
fact that the vector-current hadronic r-decay spectral functions are related 
to the isOpjvector part of the e~'"e~-annihilation cross-section via an iso-spin 
rotation l3'Q, one can obtain information from r-decay spectra to reduce 
uncertainty of e+e^ hadronic cross-sections at enecaies below the r-mass 
The inclusion of the r data, first performed in Ref.113, unfortunately lead to a 
marginal improvement only 

^"had(^^l) = 0.02810 ± 0.00062 (Alemany, Davier, Hocker 98) 

The large uncertainty of about 2.3% is due to large systematic errors of the 
total hadronic cross-sectiDiLnKjasurements. To overcome this disturbing fact 
more recent estimations MOEjO'Ej rely more on pQCD than on data. Since 
hadronic r-decays seem to be described very well by pQCD one assumes that 
perturbation theory describes the quantity Rj{s) down to Mr ~ 1.8 GeV with 
good accuracy as well, i.e., Ecut ~ 1-8 GeV is now a popular choice. The low 
energy tail below 1.8 GeV, the narrow resonances {uj, (p, J/-0's and T's) and 
the charm resonance region from 3.7 to 5.0 ate-included separately, still using 
experimental data. While the authors of Ref.ll3 use the published charm data 
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Figure 2: Different estimates of a^^ X 10 



directly, obtaining 



0.02778 ±0.00026 



(Davier, Hooker 98a) 



the authors of Refs. Eilll3 utihze the charm data only after re-normalizing in- 
dividual data sets by factors accounting for the mismatch in normalization: 
< i?!j^*''(s)/i?P*^^^(s) > averaged over small intervals justJpelow and just 
above the charmonium resonance domain. The result here isE£l 



A"had(^'4) = 0.02777 ± 0.00017 (Kiihn, Steinhauser 98) 



Another approach is based on techniques which admit su n pressing contri 
tions from data in problematic regions by a contour trickoO. This yields 



A"had(^l) = 0.02763 ± 0.00016 (Davier, Hocker 98b) 



and a~^{Mz) = 128.933 ± 0.021. 

In most cases more or less obvious theoretical uncertainties typical in hadron 
physics are not included because they are hard to be specified. The recent 
reevaluations are collected together with older results in Fig. 1. 
The second quantity of interest is the hadronic contribution to = ^ 
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which is determined from the same data by the dispersion integral 

"'^^i^j [ J — 7^ — +y — 7 — ) 

^ cut 

The kernel K{s) is a known smooth bounded function growing from 0.63... at 
threshold to 1 at oo. Note the extra 1/s-enhancement of contributions from 
low energies in Eq. as compared to Eq. ^ . Our updated evaluation reads 

ajja^ = (696.7 ± 11.9) x 10"^" (Eidelman, Jegerlehner 98) 

For this observable the use of t spectral-functions in additkui to e+e~ cross- 
sections leajis to a dramatic reduction of the uncertaintyJij. Extended use 
of pQCD i3 and application of QCD sum rule techniques t3 leads to further 
improvements 

^had ^ ^^Qj^ ^ ^ -^Q-io (Alemany, Davier, Hocker 98) 
a^'"'^ = (695.1 ± 7.5) x 10"^° (Davier, Hocker 98a) 
aji'^'i = (692.4 ± 6.2) x 10"^° (Davier, Hocker 98b) 

Again recent reevaluations are shown together with older results in Fig. 2. For 
a theory update of properties of we refer to Ref. ES. 



2 Comments on the improvements 

1) A substantial reduction of uncertainties of a^'^'^ is possible by using r-decay 
spectra to reduce the uncertainty of e+e" hadronic cross-sections at energies 
below the r-massEfEI. Conservation of the iso-vector current (CVC) allows us 
to relate non-strange r-decay data to e"'"e~-annihilation data by an iso-spin 
rotation: 

where and are related hadronic states. The e+e^ cross-section is then 
given by 

ai+l-^x'> = ^yi,x~ , V^<Mr (8) 

The r spectral function vi is obtained from the normalized invariant mass- 
squared distribution p {X~) = {1/Nx- dNx- /ds) of the r vector channel 
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X Vr by 



A 



B{t- 



B{t- 



where 



Pix-) 



(9) 



A 



Mi 



6\Vud\Hl + Skw) 



with \Vud\ = 0.9752 ± 0.0007 the CKlVLmixing matrix element and S = 0.0194 
the electroweak radiative corrections 113. With the precision of the validity of 
CVC, this allows to improve the 1=1 part of the e+e^ cross-section which by 
itself is not a directly measurable quantity. It mainly improves the knowledge 
of the TT+TT^ channel (p-resonance contribution) which is dominating in a^'^'^ 
(72%). Problems of combining data are not too serious for this particular 
observable/mode which is pure 1 = 1 and hence does not require to separate 
data into iso-spin components. The precision is then only a matter of how 
well the correspondinE_exclusive channels can be measured. With the error 
estimates discussed in Ej it was possible to reduced the error of considerably. 

In any case the following points should be kept in mind: 

(i) the usual problema-pf matching an averaging data has to be dealt with. 
For details we refer toErll3. 

(ii) In general the use of r data requires a splitting of the e+e^ data into the 
/ = and 1=1 parts. There is no precise method known how to do this, 
except for counting even and odd numbers of pions. In any case this introduces 
an additional systematic error which is hard to estimate in a clean manner. 

(iii) The previous point as well as the CVC assumption rely on estimates of 
the size of iso-spin breaking effects. This problem is absent when using e~^e~ 
data solely. Iso-spin breaking effects originate from the difference m„ — nid in 
the masses of the u and d quarks and from electromagnetic interaction, i.e., 
an imperfect treatment of QED corrections of the various exclusive channels, 
electromagnetic contributions to the hadron masses and widths etc. 

In a recent paperQ Eidelman and Ivanchenko estimate the validity of CVC by 
a detailed comparison of r channels. 



X- 



X- = tt^tt", 7r"37r°, (37r)"7r° 



with corresponding e+e channels. The results are given in the Table 1. 
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Table 1: Branching Ratios of r X v-r in percent. 



Hadronic 


World Average 


CVC 


State X 


1996 


Predictions 




25.24 ±0.16 


24.74 ±0.79 


7r-37r° 


1.14±0.14 


1.07±0.10 


(3^)-7rO 


4.25 ±0.09 


4.36 ±0.55 


Total 


30.63 ±0.23 


30.17 ± 1.00 



The difference between the CVC prediction and the measured branching ratio 
of T to 2-7r and 4-7r states is consistent with zero within the errors (-0.46 ± 
1.00)% or -1.5% relatively. There is no reason not to take this 1.5% into ac- 
count while calculating R^^^{s) from the r spectra. 

2) Substantial progress in^QCD calculations, which now includes quark mass 
effects up to three-loops EJ, allows us to apply pQCD in regions where mass 
effects are important. Many authors now assume pQCD to be valid down to 
Mt, and hence that 

CTtot(e^e^ hadrons) ^ crtot(e^e^ — > quarks and gluons) 

The assumption seems to be supported by 

• the apparent applicability of pQCD to t physics. In fact the running of 
as(MT-) as{Mz) from the r mass up to LEP energies agrees well with 
the LEP value. The estimated uncertainty may be debated, however. 

• non-perturbative (NP) effects if parametrized as prescribed by the opera^ 
tor product expansioiilDPE) of the electromagnetic current correlatorEil 
are seemingly small E3E3. 

3) Let us consider possible NP effects separately. A way to parametrize NP 
effects at sufficiently large energies and away from resonances is provided by 
the OPE aupJied to Eq. (0). Due to non-vanishing gluon and light quark 
condensatesL^ one finds the leading power corrections 

nr (g^) = ^ E Ql^^^ ■ 

q—u,d,s 



< ^GG > 
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a f 11 3, \ ,\ < mgnq > 



27 V3 486 3 "^^y / Q 

q'—u^d,s 



where a = as{^J,'^)/TT and Iq^ = \n{Q^/fi^). < ^GG > and < mqqq > are 
the scale-invariantly defined condensates. The terms beyond leading ord^r iii 
Us have been calculated from the results which were obtained in Refs. l3, cJ 
andcj. Sum rule estimates of the condensates yield typically < —GG >^ 

(0.389 GeV)4 , < rriqqq > (0.098 GeV)* for q = u,d and [< m,qq >- 

— (0.218 GeV)^ for q = s] with rather large uncertainties (see Fig. 4 below). 

Note that the expansion Eq. ( p^ ) is just a parametrization of the high energy 
tail of NP effects associated with the existence of non- vanishing condensates. 
There are other kind of NP phenomena like bound states, resonances, instan- 
tons and what else. The dilemma with Eq. (|lO|) is that it works only for 
large enough and it has been successfully applied in heavy quark physics. It 
fails do describe NP physics at lower , once it starts to be numerically rele- 
vant pQCD starts to fail because of the growth of the strong coupling constant 
(see Fig. 4 below). 



4) Some improvements are based on assuming global duality and using sum rules 
(SR). Starting point is the OPE 

< h\J^ix) MO)\h' >^Y.G^,^{x)H, 

i 

where Hi —< h\Oi{0)\h' > is a matrix element of an operator Oi between 
hadronic states h,h', and G^^i{x) are the Wilson coefficients. The global 
duality assumption asserts that all non-perturbative physics resides in the 
HiS and C^„i{x) is given by pQCD. Phenomenological tests infer that this 
works at the 10% level. Much better tests are needed to allow for precise 
confirmation as required for application in precision calculations. 
I would like to point out that the following schemes have no justification: 

* '-^had ^pQCD ^ '^resonances 

• local duality: i.e., duality applied to energy intervals (resonance regions) 
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Figure 3: Integral contour. 



Another warning I would like to make concerns pitfalls in the use of dispersion 
relations. Often one encounters arguments of the following type: consider a 
function, An'(s) = n'(s) — 11' (0) say, which is an analytic function order by 
order in perturbation theory. Analyticity then infers that the contour integral 
along a path shown in Fig. 3 vanishes. Therefore 



For a large enough circle we can apply pQCD on the right hand side and thus 
obtain the integral of our interest, which exhibits the non-perturbative physics. 
What is usually forgotten is that the uncertainty is of the order of 5 = 27ri?e 
with e being the small error expected from the truncation of the perturbative 
series. 6 easily can turn out to be large (due to R large) such that we are not 
able to make a safe estimate for the wanted integral. Since analyticity is true 
order by order in perturbation theory, we precisely reproduce the perturbative 
answer for the left hand side if we use perturbation theory on the right hand 
side. The additional use of the OPE to include the condensates, in my opinion, 
does not make the estimate much more convincing. While analyticity is a very 
powerful theoretical concept it is difficult to bo applied in numerical problems, 
because, small perturbations in one place typically cause large variations at 
remoter locations. 
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Figure 4; Adler function: data versus theory (see Ref.^*). 



3 My estimate of Aa{j j^(-M|) 



Here I propose a method for estimating Aa[j^'[j(M^) wjhich is free of assump- 
tions which are not under control. In a recent paper we investigated the 
vahdity of pQCD by means of the Adler function in the Euclidean region 
where pQCD is supposed to work best (far from resonances and thresholds). 
Results obtained from the evaluation of the dispersion integral are very well 
described by pQCD down to about 2.5 GeV. In this analysis it turned out to 
be crucial to take into account the exact mass dependence up to three loops. 
For the three-loop case this was obtained by using the series expansions pre- 
sented ino together with Pade improvement, which from a practical point of 
view yields the exact behavior in the Euclidean region. Because the precise 
mass-dependence is important we utilize the gauge iavariant background field 
MOM scheme calculation presented recently in Ref.c3. The result is shown in 
Fig. 4. With this result at hand, we then calculate in the Euclidean region 



1) 
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and obtain, for so = (2.5 GeV)^ where Aa{fj^(-so)'^^*'' = 0.007541 ± 0.000254, 

Aa|f2d(-A^l) = 0.027737 ± 0.000254 and a"i(-Mz) 128.917 ± 0.035 

for the Euchdean effective fine structure constant a value I would use in phe- 
nomcnological applications without hesitation. The virtues of this analysis are 
obvious: 

• no problems with the physical threshold and resonances 

• pQCD is used only in the Euclidean region and not below 2.5 GeV. For 
lower scales pQCD ceases to describe properly the functional dependence 
of the Adler functioncJ (although the pQCD answer remains within error 
bands down to about 1.6 GeV). 

• no manipulation of data must be applied and we need not use global or 
local duality. That contributions oLthe type Eq. (|l^) are negligible has 
been shown long time ago in Ref. E3. More recently this was confirmed 
inta. This, however, does not proof the absence of other kind of non- 
perturbative effects. 

Remaining problems are the following: 

a) contributions to the Adler function up to three-loops all have the same 
sign and are substantial. Four- and higher-orders could still add up to non- 
negligible contribution. An error for missing higher order terms is not included. 
The scheme dependence MS versus background field MOM has been discussed 
in Ref.0. 

b) the link between space-like and time-like region is the difference 

A = Aa'^^^iMj) - A4ad(-Af|) = 0.000045 ± 0.000002 , 

which can be calculated in pQCD. It accounts for the Z7r-terms from the logs 
ln(— (7^//i^) = ln(|g^//i^|) +iTr. One may ask the question whether these terms 
should be resummed at all, i.e., included in the running coupling. Usually such 
terms tend to cancel against constant rational terms which are not included 
in the renormalization group (RG) evolution. It is worthwhile to stress here 
that the running coupling is not a true function of (or even an analytic 
function of ) but a function of the RG scale fi^ . The coupling as it appears 
in the Lagrangian in any case must be a constant, albeit a /x^-dependent one, 
if we do not want to end up in conflict with basic principles of quantum field 
theory. The effective identification of iJ? with a particular value of q^ must be 
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understood as a subtraction (reference) point. 

The above result was obtained using the gauge invariapt background-field 
MOM renormalization scheme, presented recently in Ref. EZl. In the transition 
from the MS to the MOM scheme we adapt the rescaling procedure described 
in £3, such that for large /i 

as{M^)=as{fi^) + + Oia^,) . (12) 

This means that xq is chosen such that the couplings coincide to leading and 
next -to-leading order at asymptotically large scales. Numerically we find xq ~ 
2.0144. Due to this normalization by rescaling the coefficients of the Adler- 
function remain the same in both schemes up to three-loops. In the MOM 
scheme we automatically have the correct mass dependence of full QCD, i.e., we 
have automatic decoupling and do not need decoupling by hand and matching 
conditions likcj-Ln the MS scheme. For the numerical evaluation we use the pole 
quark masses N = 1.55GeV, nn, = 4.70GeV, rrit = 173.80GeV and the 
strong interaction coupling a^^- — (Mz) — 0.120 ±0.003. For further details we 
refer to El. 



4 Improvements expected from Daphne 

For the rather dramatic progress on the theory-driven reduction of uncertain- 
ties of Aa[^^{M'^) and a^, discussed in Sec. 1, little compelling theoretical or 
experimental justification is available to remove doubts concerning the exis- 
tence of further unaccounted theoretical uncertainties. Therefore it is impor- 
tant to pursue the experimental program to reduce uncertainties of R~f{s) mea- 
surements. New low energy data are expected from BES, CMD II and DA$NE 
in the near future. On a time scale of a few years the <I>-factory Daphne is the 
most promising facility. Therefore one of the challenges for Daphne is presented 
by the measurement of the total e+e~-hadronic cross-section at y/s^2 GeV, 
which provides the information on the hadronic contributions to — and to 
a lesser extent to a{Mz)- Considerations from precision physics suggest that 
the goal should be to measure the cross-section at the level of a few permille. 
This could allow a reduction of the error on such as tii be closer to the 
contribution from virtual light-by-light scattering (see e.g.E3). To extract the 
total cross-section a certain amount of modeling is necessary, since, in contrast 
to LEP, at low energies the theoretical shape of the total cross-section cannot 
be calculated from first principles. Important cross-checks will be obtained 
by studying connections between form-factor measurements at Daphne and r- 
lepton decays (CLEO, ALEPH, OPAL, DELPHI, BABAR) through iso-spin 
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Figure 5: Multi-particle channels in e+e -annihilation at low energy. 

relations between multi-pion final states, wliic;li wc; mentioned before. 

There are two stages: at the <l?-resonance one can measure the pion form 
factor |i^7r(s)| for ^/s ^ 1 GeV via hard photon tagging 

$ ^ 7* + 7(real, hard) n^n^ + 7(real, hard) 

At a second stage, after a non-trivial Daphne upgrade, one expects to be able 
to perform an energy scan from threshold up to about 2 GcV, which allows 
for a "direct" measurement of i?.^(.s). Unfortunately the measurement of this 
inclusive quantity at low energies is not as easy as at LEP/SLD. Because of 
competing 2 body channels e+e^, /x+/i^, tt+tt^, K^K^, ■■■ , which must 
be unambiguously separated, a good particle identification is required. This 
means that the inclusive cross-section can be obtained only by summing up 
the exclusive channels. This also requires a detailed theoretical understanding 
of the individual channels at low energies (^2 GeV). In Fig. 5 we show the 
spectrum of channels which contribute the total hadronic cross-section. The 
contributions from the exclusive channels to a]','^'^ are listed in Table 2. In order 
to reduce the uncertainty from the region below 2 GeV to the 0.3% level the 
accuracy required for the individual channels is given in the last column. 



m oj m 4> 

I — H — I p 

O.bo 500'. 00 1000.00 1500.00 
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Table 2: Contributions to d,, = a,, x 10^^ from exclusive channels. 
A* M 



channel 


-had 


acc. 


channel 


-had 


acc. 




506 


0.3% 


Stt 


4 


10% 




47 


- 1% 




4 


i 




40 


i 


KsKl 


1 






24 






1.8 






14 






0.5 






5 


10% 




0.2 




2 GeV < £; < Afj/^ 


22 










Mj/^ <E< Mr 


20 










Mr < E 


<5 











Theoretical work with the aim to calculate radiative corrections at the level of 
precision as indicated in Table 2 is in progress. 

5 Summary and conclusions 

The future of precision physics depends to a large extent on the uncertainties 
of hadronic vacuum polarization effects. Present and future are different for 
flhad and Aa(t)^(M|): 

^had . uncertainty from TT+TT^. jo'^ dominates 



Table 3: Uncertainties of a^^'^ in units 10 



5a^ 


input 


5a^ 


input 


- 156 

- 60 

- 62 


e+e- data&0 

Daphne + additional r data 

theory-driven guess El 


- 102 

- 40 


e+e-, T data0 

BNL a^-experimentlii 



^°^had(-^I) • uncertainty distributed everywhere below T 

Uncertainties of e+e^i-data range from 5% to 20% at present which yields 
(5Aq ~ 0.00065 (data)Q't2l. Suppose we could improve this to a 1% measure- 
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ment in systematics and with high enough statistics, this would lead us to 
SAa ^ 0.00028 (data), whichJaasL to be confronted with the theory driven es- 
timate SAa - 0.00016 (theo)li3M There is little hope that this precision can 
ever be reached by an experiment. 
Alternative strategies: 

• try to determine Aa^^^^{M'^) by a direct measurement 

• simplest: find the Higgs and determine mn precisely; this would allow 

us to determine Aa^^^{M^) from the precise value of sin^ 9^^^ measured 
at LEP/SLD. 

• For what concerns theory, which is essentially pQCD, we need better 
estimates of uncertainties (scheme dependence, higher orders_etc. which 
are a problem at Mr), exterided studies in MOM schemes cll, analytic 
extensions of the MS scheme l2I etc. are needed. 

In the meantime I propose to use in phenomenological applications an estimate 
of the kind I presented in Sec. 3. 
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